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[i]  The  7  years  (2003-2009)  of  output  from  a  regional  version  of  the  Hybrid  Coordinate 
Ocean  Model  (HYCOM)  nested  in  global  HYCOM  are  used  to  study  the  seasonal  and 
interannual  variability  of  the  salinity  in  the  northern  Gulf  of  California  (NGOC).  Previous 
studies  illustrate  that  the  NGOC  is  characterized  by  an  annual  evaporation  of  ~0.9  m/yr.  This 
evaporation  generates  high-sea-surface-salinity  (SSS)  water,  which  reaches  a  maxima  in 
the  NGOC  (>37)  and  decreases  to  ~35  toward  the  entrance  of  the  Gulf  of  California.  The 
NGOC  surface  water  is  interannually  modulated  by  fluctuations  in  evaporation  and  by 
fluctuations  in  the  low-salinity  water  transported  into  the  region  by  poleward  eastern 
boundary  currents.  The  fluctuations  in  the  transport  of  low-salinity  water  are  linked  to  the 
arrival  of  equatorially  originated  coastally  trapped  waves.  The  crucial  role  of  the  transport  of 
low-salinity  water  for  the  interannual  variability  of  SSS  is  illustrated  by  the  2006  and 
2008  fall  seasons,  which  include  the  lowest  simulated  salinity  of  the  period  2003-2009.  The 
lowest  salinity  in  2006  and  2008  cannot  be  explained  solely  by  evaporation  because  2006 
was  characterized  by  the  largest  evaporation  of  the  period  2003-2009.  However,  the 
presence  of  the  lowest-salinity  environment  can  be  attributed  to  the  evaporation  in 
conjunction  with  the  largest  upper  ocean  transport  of  low-salinity  water  carried  to  the  NGOC 
by  the  2006  and  2008  coastally  trapped  waves  intensified  poleward  eastern  boundary  currents. 
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1.  Introduction 

[2]  The  northern  Gulf  of  California  (NGOC)  is  a  shallow 
and  semienclosed  basin  characterized  by  high-evaporation 
(0.9-0.95  m/yr)  and  high-sea-surface-salinity  (SSS)  waters 
(>37.0)  [Bray,  1988;Tavi«  and Organista,  1988;Lflv/«  etal, 
1998].  The  high  solar  insolation,  geographical  isolation  of  the 
NGOC  (which  is  surrounded  by  deserts  and  mountains)  and 
the  arrival  of  several  cold  fronts  that  transport  cold,  dry  desert 
air  into  the  region,  all  combine  to  produce  high  evaporation 
[Bray,  1988;  Lavin  and  Organista,  1988;  Badan-Dangon 
etal,  1991;  Reyes  and  Lavin,  1997;  Lopez,  1997]. 

[3]  Based  on  hydrographic  observations,  Alvarez -Borrego 
and Schwartzlose  [\919'\,Bray  [1988],  andLavin  etal.  [1995, 
1998]  hypothesized  that  the  high  evaporation  in  the  NGOC 
would  produce  high-SSS  water.  In  particular,  the  data 
reported  by  Lavin  et  al.  [1998]  shows  the  highest  salinity  in 
the  shallow  upper  NGOC.  In  the  same  line  of  research,  Lopez 
[1997]  used  a  numerical  ocean  model  to  study  the  NGOC 
winter  circulation  and  finds  that  as  a  result  of  cooling  and 
evaporation,  water  in  the  shallow  NGOC  increases  its  den¬ 
sity.  Emphasizing  the  importance  of  the  oceanic  remote 
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forcing,  Ripa  [1997]  reported  that  most  of  the  dynamics  and 
thermodynamics  of  the  Gulf  of  California  (GOC)  are  con¬ 
trolled  by  the  Pacific  Ocean,  which  excites  a  baroclinic 
Kelvin  wave  at  the  entrance  of  the  GOC.  Based  on  the  Ripa’s 
[1997]  Kelvin  wave  scenario,  Berdn-Vera  and  Ripa  [2002] 
used  a  linear  one-dimensional  nondiffusive  inhomogeneous 
two-layer  model  to  study  the  salt  balance  in  the  GOC  and 
found  that  the  seasonal  balance  of  the  salinity  is  largely 
controlled  by  the  Pacific  Ocean. 

[4]  In  addition  to  the  several  questions  answered  by  these 
studies,  some  valid  questions  have  not  yet  been  answered. 
Is  the  NGOC  high  evaporation  the  essential  generator  and 
modulator  of  the  NGOC  high-salinity  water?  What  is  the  role 
of  the  Costa  Rica  Coastal  Current  (CRCC)  and  the  Mexican 
Coastal  Current  (MCC)  on  the  seasonal  evolution  of  the 
NGOC  upper  ocean  water  and  upper  ocean  circulation?  What 
is  the  role  of  equatorially  forced  coastally  trapped  waves 
(CTWs)  on  the  interannual  variability  of  the  NGOC  surface 
water? 

[5]  This  work  builds  on  the  previous  studies  by  doc¬ 
umenting  the  NGOC  seasonal  evolution  of  SSS  water,  which 
has  a  strong  interannual  variability,  and  is  modulated  by  the 
NGOC  high  evaporation,  and  by  the  CRCC  and  the  MCC.  To 
explain  the  origin  of  the  seasonal  and  interannual  salinity 
variability  of  the  NGOC  surface  water,  two  hypotheses  are 
presented  and  discussed.  The  first  hypothesis  is  that  evapo¬ 
ration  is  the  primary  generator  and  modulator.  The  second 
hypothesis  builds  on  Ripa’s  [1997]  model  of  the  Pacific 
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Table  1.  Summary  of  the  Main  Features  of  the  Simulations  Used  in  This  Study’* 


Experiment  Domain 

Boundary  Conditions 
From  Experiment 

Horizontal 

Resolution 

Period 

(years) 

Winds 

Wind  Remote 

Frequency  Forcing  CRCC  MCC 

Annual  ^ 
CTW 

Interannual 

CTW 

1 

Global 

N/A 

1/25° 

1-10 

ERA40/QuikSCAT 

Monthly 

N/A 

Yes 

Yes 

Yes 

No 

2 

Global 

N/A 

1/25° 

2003-2009  NOGAPS/QuikSCAT 

6  h 

N/A 

Yes 

Yes 

Yes 

Yes 

3 

Nested  GOC 

1 

1/25° 

1-10 

ERA40/QuikSCAT 

Monthly 

Yes 

Yes 

Yes 

Yes 

No 

4 

Nested  GOC 

1 

1/25° 

1-10 

No  Wind 

Monthly 

Yes 

Yes 

No 

Yes 

No 

5 

Nested  GOC 

N/A 

1/25° 

1-10 

ERA40/QuikSCAT 

Monthly 

No 

No 

Yes 

No 

No 

6 

Nested  GOC 

2 

1/25° 

2003-2009  NOGAPS/QuikSCAT 

6  h 

Yes 

Yes 

Yes 

Yes 

Yes 

7 

Nested  GOC 

2 

1/25° 

2003-2009 

No  Wind 

6  h 

Yes 

Yes 

No 

Yes 

Yes 

8 

Nested  GOC 

N/A 

1/25° 

2003-2009  NOGAPS/QuikSCAT 

6  h 

No 

No 

Yes 

No 

No 

“All  the  experiments  included  in  this  study  contain  32  layers. 


Ocean  as  the  main  forcing  of  the  GOC.  Thus,  together  these 
two  hypotheses  suggest  that  the  evaporation  in  conjunction 
with  the  low-salinity  water  (~35)  transported  by  the  CROC, 
the  MCC,  and  the  CTWs  are  responsible  for  the  seasonal  and 
interannual  variability  of  the  SSS.  Our  approach  is  based  on 
the  analysis  of  the  results  of  three  climatological  and  three 
interannual  forced  regional  GOC  Hybrid  Coordinate  Ocean 
Model  (HYCOM)  simulations,  which  are  nested  in  global 
HYCOM.  The  nested  approach  allows  the  isolation  of  the 
effects  of  both  local  and  remote  forcings  and  the  time  period 
of  the  interannual  simulations  (2003-2009)  allows  the  eval¬ 
uation  of  the  effects  of  the  El  Nino-La  Nina  cycle  on  the 
NGOC. 

2.  The  Model 

[6]  The  Hybrid  Coordinate  Ocean  Model  is  isopycnal  in 
the  open  stratified  ocean,  terrain  following  in  shallow  coastal 
regions,  and  maintains  z  levels  in  the  mixed  layer  and 
unstratified  regions.  This  generalized  vertical  coordinate 
approach  is  dynamic  in  space  and  time  via  the  layered  con¬ 
tinuity  equation  that  allows  a  smooth  dynamical  transition 
between  the  coordinate  types.  HYCOM  [Bleck,  2002]  was 
developed  from  the  Miami  Isopycnic  Coordinate  Ocean 
Model  using  the  theoretical  foundation  for  implementing 
a  hybrid  coordinate  system  [Bleck  and  Benjamin,  1993]. 
Application  of  HYCOM  to  the  GOC  has  been  discussed  by 
Lopez  et  al.  [2005]  and  Zamudio  et  al.  [2008,  2010]. 

[?]  Two  different  eddy-resolving  (1/25°  equatorial  reso¬ 
lution)  HYCOM  configurations  are  used.  The  global  model 
(which  is  the  provider  of  boundary  conditions  for  the  nested 
GOC  model)  is  on  a  Mercator  projection  from  78.64°S  to 
47°N  and  north  of  this  it  employs  an  Arctic  dipole  patch 
where  the  poles  are  shifted  over  land  to  avoid  a  singularity  at 
the  North  Pole.  The  nested  GOC  model  is  on  a  Mercator 
projection  and  extends  from  118°W  to  97° W  and  14°N  to 
32°N.  The  global  model  was  initialized  using  temperature 
and  salinity  from  the  Generalized  Digital  Environmental 
Model  (GDEM)  Version  3.0  climatology  [Carnes,  2009]. 
Next,  it  was  spun-up  for  10  years  using  a  monthly  climatol¬ 
ogy  based  on  the  1979-2002  1.125°  European  Centre  for 
Medium-Range  Weather  Forecasts  (ECMWF)  40  year 
reanalysis  (ERA40)  [Kallberg  et  al,  2004]  using  the  bulk 
formulation  of  Kara  et  al.  [2005a]  for  converting  10  m  winds 
to  wind  stress  and  with  the  wind  speeds  corrected  using  a 
regression  to  a  monthly  climatology  from  QuikSCAT  scat- 
terometer  data  [Kara  et  al,  2009],  i.e.,  ERA40-QuikSCAT, 
as  previously  implemented  by  Metzger  et  al.  [2010].  In 
addition,  the  model  is  driven  by  daily  averaged  surface  air 


temperature,  surface  atmospheric  humidity,  surface  short¬ 
wave  and  long-wave  heat  fluxes,  and  precipitation.  Surface 
latent  and  sensible  heat  fluxes  and  evaporation  are  calculated 
at  each  model  time  step  using  the  model’s  top  layer  temper¬ 
ature  [Kara  et  al.,  2005a].  Furthermore,  the  model  includes 
monthly  rivers  [Barron  and  Smedstad,  2002]  and  turbidity 
forcing  [Kara  et  al.,  2005b],  realistic  bottom  topography  and 
coastline  geometry,  which  are  based  on  a  modified  version  of 
the  NRL  DBDB2  topography  (http://www7320.nrlssc.navy. 
mil/DBDB2_WWW).  The  model  uses  the  5  m  isobath  as  a 
land-sea  boundary,  includes  32  vertical  coordinate  layers, 
allows  isopycnals  to  intersect  sloping  topography  by  allowing 
zero  thickness  layers,  does  not  include  ocean  data  assimila¬ 
tion,  and  relaxes  the  sea  surface  salinity  toward  the  Polar 
Science  Center  Hydrographic  Climatology  [Steele  et  al., 
2001].  This  is  required  to  prevent  long-term  SSS  drift 
due  to  inadequate  precipitation  minus  evaporation  forcing. 
After  spin-up,  the  model  used  the  Fleet  Numerical  Meteo¬ 
rology  and  Oceanography  Center  (FNMOC)  3  hourly  0.5° 
Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  forcing  [Rosmond  etal,  2002]  spanning  the  time 
frame  2003-2009.  For  the  wind  forcing,  the  long-term 
(annual)  mean  from  NOGAPS  was  replaced  by  ERA40  mean 
for  consistency  on  the  larger  scale  as  discussed  by  Metzger 
etal.  [2010]. 

[s]  The  nested  GOC  H Y COM  was  similarly  initialized  and 
forced  as  the  global  model  in  both  climatological  and  inter¬ 
annual  modes.  However,  lateral  boundary  conditions  for  the 
nested  model  were  provided  by  global  HYCOM.  In  addition, 
because  the  Colorado  river  in  the  NGOC  region  has  been 
dammed  [e.g.,  Carbajal  et  ah,  1997;  Lavin  and  Sanchez, 
1999],  the  nested  model  does  not  include  any  river  input 
and  the  SSS  is  not  relaxed  to  any  climatology.  The  compar¬ 
ison  of  model  results  and  observations  of  section  3 . 1  validates 
this  approach.  All  the  simulations  discussed  in  the  present 
study  are  included  in  Table  1.  Furthermore,  section  3.2.1 
describes  the  application  of  the  different  nested  simula¬ 
tions  to  separate  the  circulation  induced  by  local  and  remote 
forcings. 

3.  Results  and  Discussion 

3.1.  Annual  Salinity 

[9]  The  long-term  mean  SSS  of  the  GOC  is  characterized 
by  a  minimum  of  ~35  at  the  entrance  of  the  GOC  (Figure  la). 
That  increases  toward  the  NGOC  and  reaches  the  maximum 
of  ~36.6  at  the  upper  GOC  and  incorporates  the  inverse 
estuarine  features  of  the  upper  NGOC  (maximum  salinity  at 
the  head  of  the  NGOC),  which  were  previously  documented 
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Figure  1.  (a)  The  5  year  mean  sea  surface  salinity  (color  contours)  as  determined  by  a  climatologically 
forced  nested  GOC  HYCOM  simulation  (experiment  3  in  Table  1).  (b)  Sea  surface  salinity  for  23  June 
to  1  July  1996.  Adapted  from  Lavin  et  al.  [1998],  with  permission  from  Elsevier,  (c)  Sea  surface  salinity 
(color  contours)  and  upper  ocean  (0-100  m)  currents  (arrow  vectors)  for  23  June  2006  as  determined  by  an 
interannually  forced  nested  GOC  HYCOM  simulation  (experiment  6  in  Table  1).  The  upper  ocean  transport 
of  Figure  6c  was  calculated  along  the  white  zonal  line  located  at  the  entrance  of  the  GOC.  (d)  Meridional 
currents  (color  contours  in  cm/s)  and  (e)  salinity  for  23  June  2006  along  the  white  zonal  line  at  the  entrance 
of  the  GOC  in  Figure  Ic. 


by  Lavin  et  al.  [1998]  (Figure  lb).  Although  the  model 
simulates  the  observed  inverse  estuarine  features  of  the  upper 
NGOC,  the  comparison  of  the  SSS  ofFigures  la  and  lb  might 
suggest  that  the  model  underestimates  the  magnitude  of  the 


maximum  observed  SSS  of~37.0.  However,  it  is  important  to 
keep  in  mind  that  Figure  1  a  depicts  a  5  year  mean  of  simulated 
SSS  whereas  Figure  lb  depicts  the  instantaneous  observed 
SSS  during  a  few  days  in  early  summer  of  1996.  So,  to  pro- 
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vide  some  insight  on  the  validation  of  the  model  results, 
Figure  Ic  shows  a  simulated  instantaneous  SSS  field  during 
early  summer  of 2008  that  features  a  maximum  SSS  of~37.1 
and  this  compares  better  with  the  maximum  observed  SSS  of 
~37.0.  Even  though  the  simulations  of  Figure  1  do  not  include 
tides,  and  consequently  they  do  not  include  tidal  mixing 
(which  has  been  proposed  as  an  important  mechanism 
affecting  the  distribution  of  the  salinity  in  the  NGOC  [e.g., 
Argote  et  al,  1995;  Lavin  et  al,  1998]),  they  simulate 
the  general  characteristics  of  the  GOC  SSS  reported  in  the 
literature. 

[10]  The  instantaneous  SSS  field  displays  low-salinity 
water  transported  to  the  north  by  the  poleward  eastern 
boundary  currents  (PEBC)  along  the  west  coast  of  mainland 
Mexico  (Figure  Ic).  During  this  day  the  PEBC  extend  from 
the  surface  to  ~50  m  depth  and  advect  water  with  salinity 
values  <35  (Figures  Id  and  le).  The  low-salinity  tongue 
extends  from  the  entrance  of  the  GOC  (with  salinity  <3  5)  to 
the  NGOC  (with  salinity  of  >37.0)  and  during  its  northward 
advection  is  modified  (increasing  its  salinity)  by  the  water  of 
the  GOC.  In  addition  to  a  strong  latitudinal  gradient  (which 
increases  northward,  reaches  the  maximum  in  the  NGOC,  and 
it  is  included  in  both  the  instantaneous  and  the  long-term 
field),  the  instantaneous  SSS  field  is  characterized  by  a  lon¬ 
gitudinal  gradient  that  increases  westward  and  reaches  the 
maximum  along  the  eastern  coast  of  the  Baja  California 
Peninsula  (Figure  Ic)  as  measured  by  Lavin  et  al.  [1998]  and 
Castro  et  al.  [2000].  This  salinity  maximum  is  not  limited  to 
the  sea  surface,  it  extends  to  ~125  m  depth  and  it  is  advected 
southward  by  a  narrow  western  boundary  equatoward  current 
(Figures  Id  and  le). 

3.2.  Monthly  Salinity  and  Circulation 

3.2.1.  Separation  of  Circulation  Induced  by  Local 

and  Remote  Forcings 

[11]  Considering  that  monthly  maps  for  the  six  different 
GOC  simulations  of  Table  1  would  require  six  different 
12  panel  figures,  the  separation  of  the  circulation  induced  by 
local  and  remote  forcings  is  presented  and  discussed  using  the 
month  of  July  as  a  representative  example.  Figure  3,  which 
includes  12  months  of  the  NGOC  circulation,  is  shown  and 
discussed  in  section  3.2.2. 

[12]  Figure  2  shows  upper  ocean  PEBC  from  Cabo  Cor- 
rientes  (in  the  southeast  comer)  to  the  NGOC.  These  currents 
reach  the  maximum  (secondary  maximum)  northward  pene¬ 
tration  into  the  NGOC  during  July  (December-January) 
producing  the  primary  (secondary)  annual  maximum  in  the 
northward  transport  [see  Zamudio  et  al,  2008,  Figure  5].  The 
currents  of  Figure  2  are  composed  of  (1)  the  CRCC  that 
originates  from  a  branch  of  the  North  Equatorial  Counter 
Current  that  turns  northward  near  the  western  coast  of  the 
Americas  [Wyrtki,  1966;  Zamudio  et  al,  2001;  Kessler, 
2006],  (2)  the  currents  associated  with  the  propagation  of 
low-frequency  CTWs  [Ripa,  1997;  Strub  and  James,  2002a, 
2002b],  and  (3)  the  MCC  that  originates  as  the  oceanic 
response  to  the  local  wind  stress  curl  [Zamudio  et  al.,  2007; 
Godinez  et  al,  2010].  Isolating  the  physical  processes  that 
contribute  to  the  variability  of  the  NGOC  SSS  is  one  of  the 
main  objectives  of  this  study.  Because  the  PEBC  are  pro¬ 
posed  as  the  main  transporters  of  low-salinity  water  from  the 
Pacific  to  the  NGOC,  accuracy  in  the  simulation  of  the  PEBC 
is  essential  to  simulate  the  GOC  circulation,  including  the 


NGOC  seasonal  reversal  of  the  circulation,  which  was  pre¬ 
viously  observed  by  Lavin  et  al.  [1997]  and  modeled  by  Beier 
and  Ripa  [1999]  and  Marinone  [2003]  who  used  analytical 
wind  and  CTWs  to  force  regional  models.  Flere,  we  build 
upon  these  valuable  model  studies  hypothesizing  that  the 
three  components  of  the  PEBC  contribute  to  the  generation  of 
the  NGOC  seasonal  reversing  circulation.  To  test  this 
hypothesis,  the  effects  of  the  CRCC  and  the  annual  CTW  on 
the  generation  of  the  NGOC  seasonal  reversing  circulation 
are  isolated  by  forcing  the  regional  model  with  oceanic 
remote  forcing  only.  Next,  the  effects  of  the  MCC  on  the 
formation  of  the  NGOC  seasonal  reversing  circulation  are 
isolated  by  forcing  the  regional  model  with  local  forcings 
only.  Finally  a  third  simulation  is  forced  with  both  local  and 
remote  forcings.  The  different  currents  and  CTWs  included  in 
these  simulations  are  summarized  in  Table  1. 

[13]  There  are  several  salient  features  of  the  simulations  in 
Figure  2.  All  the  simulations  consist  of  upper  ocean  PEBC 
from  Cabo  Corrientes  to  the  NGOC,  a  cyclonic  circulation  in 
the  NGOC,  and  a  series  of  eddies  along  the  axis  of  the  GOC. 
The  two  simulations  forced  with  local  plus  remote  forcing 
(Figures  2a  and  2b)  include  the  strongest  and  most  continuous 
currents  of  all  the  simulations.  This  suggests  positive  cou¬ 
pling  (strengthening)  between  the  currents  generated  by  local 
and  remote  forcings  since  the  simulations  forced  with  remote 
forcing  only  (Figures  2c  and  2d)  and  local  forcing  only 
(Figures  2e  and  2f)  generate  by  themselves  weaker  summer 
GOC  PEBC  and  weaker  NGOC  cyclonic  circulation  than  the 
currents  in  Figures  2a  and  2b.  In  addition  to  the  PEBC,  the 
simulations  in  Figures  2a-2d  incorporate  equatorward  west¬ 
ern  boundary  currents  along  the  western  coast  of  the  GOC, 
from  the  NGOC  to  the  southern  tip  of  the  Baja  California 
Peninsula.  These  currents  can  mainly  be  attributed  to  the 
remote  forcing  because  the  equatorward  western  boundary 
currents  are  not  evident  in  the  simulations  lacking  remote 
forcing  (Figures  2e  and  2f).  Furthermore,  there  are  important 
differences  in  the  magnitude  of  the  currents  modeled  in  the 
climatological  and  interannual  simulations  (Figure  2a  versus 
Figure  2b,  Figure  2c  versus  Figure  2d,  and  Figure  2e  versus 
Figure  2f),  which  can  be  attributed  to  the  effects  of  wind 
events  and  the  interarmual  variability  of  the  remote  forcing.  In 
short,  both  the  local  and  the  remote  forcings  act  individually 
or  together  to  induce  upper  ocean  PEBC  that  penetrate 
into  the  NGOC  and  generate  the  observed  summer  cyclonic 
circulation. 

3.2.2.  NGOC  Monthly  Salinity  and  Circulation 

[14]  The  monthly  upper  ocean  currents  of  Figure  3  reveal  a 
strong  seasonal  fluctuation  in  the  circulation  of  the  NGOC. 
That  is  distinguished  by  an  overall  seasonal  reversal  of  the 
circulation.  From  November  to  April  anticyclonic  currents 
dominate  the  circulation  in  the  NGOC.  This  anticyclonic 
circulation  evolves  from  occupying  the  complete  NGOC, 
including  the  circulation’s  center  close  to  113.6°W,  30.3°N 
and  maximum  current  speed  of  ~15  cm/s  during  November 
(Figure  3k),  to  a  more  concentrated  circulation  with  the  center 
close  to  113.8°W,  30°N  and  maximum  current’s  speed  of 
~1 1  cm/s  in  April  (Figure  3d).  Furthermore,  from  November 
to  April  there  is  conclusive  evidence  of  the  effect  of  the  PEBC 
on  neither  the  anticyclonic  circulation  nor  the  SSS  (Figure  4). 
That  is  an  expected  result,  since  the  NGOC  transport  induced 
by  the  PEBC  reaches  the  annual  minimum  during  November 
[see  Zamudio  et  al,  2008,  Figure  5]  and  consequently  the 
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Figure  2.  July  mean  upper  ocean  (0-100  m)  currents  as  determined  by  six  different  nested  GOC  HYCOM 
simulations.  The  color  contours  (in  cm/s)  represent  the  magnitude  of  the  currents  and  the  arrow  vectors  rep¬ 
resent  the  direction.  The  experiment  number  (described  in  Table  1)  is  indicated. 


simulated  SSS  minimum  is  >35.8  from  November  to  March 
(Figure  4).  In  addition,  from  November  to  March  the  NGOC 
is  characterized  by  a  general  southeastward  directed  wind 
(Figure  4),  which  opposes  the  propagation  of  the  PEBC 
weakening  the  northward  penetration.  Nevertheless,  note  the 


penetration  of  the  PEBC  to  the  NGOC  during  January  that 
is  consistent  with  the  December-January  secondary  maxi¬ 
mum  of  the  northward  transport  [see  Zamudio  et  ah,  2008, 
Figure  5].  In  April  the  anticyclonic  circulation  starts  to 
weaken  due  to  the  variation  in  the  wind’s  direction  that 
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Figure  3.  Monthly  mean  upper  ocean  (0-100  m)  currents  as  determined  by  a  climatologically  forced 
nested  GOC  HYCOM  simulation  (experiment  3  in  Table  1).  The  color  contours  (in  cm/s)  represent  the 
magnitude  of  the  currents  and  the  arrow  vectors  represent  the  direction.  White  polygons  are  model  land. 


changes  from  southeastward  directed  to  more  eastward 
directed  (Figures  3d  and  4d).  By  May  the  PEBC  starts  to 
penetrate  into  the  NGOC  diminishing  the  anticyclonic  cir¬ 
culation  (Figure  3e).  Thus,  the  transition  month  of  May  is 
characterized  by  a  meandering  poleward  current.  From  June 
to  September  the  wind  blows  in  a  general  northeastward 
direction  (Figures  4f-4i)  and  the  NGOC  features  an  upper 
ocean  cyclonic  circulation  (Figures  3f-3i).  During  these 
summer  months  the  upper  ocean  currents  reach  the  maximum 
speed  (~40  cm/s)  for  the  year  and  the  communication  of 
the  NGOC  circulation  with  the  southern  GOC  is  more  pro¬ 
nounced  than  during  fall-winter,  as  evidenced  by  the  low- 
SSS  tongue  extending  cyclonically  from  along  the  coast  from 
June  to  September  (Figure  4f-4i).  During  the  course  of 
September  the  cyclonic  circulation  starts  to  weaken  and  by 


the  transition  month  of  October  the  NGOC  is  characterized  by 
a  meandering  equatorward  circulation  (Figure  3j). 

[is]  Overall,  Figures  3  and  4  show  that  the  influence  of  the 
Pacific  Ocean  and  the  southern  GOC  on  the  NGOC  upper 
ocean  circulation  and  SSS  is  stronger  during  spring-summer 
than  during  fall-winter.  Furthermore,  independent  of  the 
monthly  SSS  variability  the  maximum  SSS  is  always  located 
at  the  head  of  the  upper  NGOC  and  it  reaches  the  yearly 
maximum  during  summer  as  measured  by  Lavin  etal.  [1998]. 
To  provide  some  insight  into  the  role  of  the  PEBC  as  trans¬ 
porters  of  low-salinity  water  from  the  Pacific  to  the  NGOC, 
we  calculate  the  monthly  maximum  of  the  NGOC  SSS  for  a 
simulation  forced  with  local  and  remote  forcings  (thick  line  in 
Figure  5),  and  for  a  simulation  forced  with  local  forcings  only 
(thin  line  in  Figure  5).  The  two  simulations  capture  the  annual 
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Figure  4.  Monthly  mean  sea  surface  salinity  (color  contours)  as  determined  by  a  climatologically  forced 
nested  GOC  HYCOM  simulation  (experiment  3  in  Table  1)  and  10  m  wind  stress  (arrow  vectors)  simulated 
by  NOGAPS.  White  polygons  are  model  land. 


cycle,  but  the  simulation  forced  with  local  forcing  only  pro¬ 
duces  higher  SSS  than  the  simulation  forced  with  local  plus 
remote  forcings.  A  plausible  explanation  for  this  result  is 
as  follows.  The  high  evaporation  in  the  NGOC  generates 
high-SSS  water  and  the  PEBC  modulates  the  generation 
of  the  high-SSS  water  by  transporting  upper  ocean  low- 
salinity  water  from  the  Pacific  to  the  NGOC  as  shown  in 
Figure  Ic.  The  high  evaporation  and  the  MCC  are  included  in 
the  two  simulations  of  Figure  5a.  However,  the  CTWs  and  the 
CRCC  (which  transport  low-salinity  water  from  the  coastal 
regions  of  Southern  Mexico  to  the  GOC  (Animation  1  *  and 
section  3.3))  are  only  incorporated  in  the  simulation  with 


'Animation  1  is  available  in  the  HTML. 


remote  forcings  (thick  line  in  Figure  5a).  Furthermore,  it  is 
important  to  recognize  the  dominance  of  the  local  forcing  on 
the  NGOC  SSS  monthly  maximum  regardless  of  the  influ¬ 
ence  of  the  PEBC  freshwater  transport  (Figure  5a).  Equally 
significant  are  the  effects  of  the  PEBC  freshwater  transport  on 
amplifying  the  semiannual  cycle  of  the  NGOC  SSS  monthly 
mean  (Figure  5b). 

3.3.  Interannual  Variability  of  the  Salinity 

[16]  The  NGOC  HYCOM  simulated  SSS  for  the  period 
2003-2009  shows  strong  interannual  variability  (Figure  6a). 
As  a  response  to  the  high  evaporation  in  the  NGOC 
(Figure  6b),  water  with  salinities  >36  is  generated  at  the  sea 
surface  during  the  7  years  of  this  study  (Figure  6a).  Indeed,  in 
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Figure  5.  Time  series  of  (a)  monthly  maximum  and  (b) 
monthly  mean  sea  surface  salinity  for  the  NGOC  region 
shown  in  Figure  4,  as  determined  by  a  climatological  simu¬ 
lation  forced  with  local  plus  remote  forcing  (thick  line  labeled 
as  L  +  RF  is  experiment  3  in  Table  1)  and  by  a  climatological 
simulation  forced  with  local  forcing  only  (thin  line  labeled  as 
LF  is  experiment  5  in  Table  1). 

2005  the  maximum  SSS  reaches  salinities  >36.3,  and  during 
2003  and  2007  it  reaches  salinities  >36.2.  In  contrast,  during 

2006  and  2008  the  minimum  SSS  reaches  salinities  <35.7 
(Figure  6a).  Why  do  2005  and  2007  have  the  saltiest  water  of 
the  period  2003-2009?  Why  are  2006  and  2008  characterized 
by  the  freshest  water?  Answers  to  these  questions  are  pre¬ 
sented  and  discussed  in  sections  3.3.1  and  3.3.2. 

3.3.1.  Evaporation  as  the  Essential  Generator 
and  Modulator  of  the  NGOC  Salinity 
[17]  Assuming  that  the  interannual  variability  in  the  gen¬ 
eration  of  the  NGOC  high-salinity  water  is  predominantly 
due  to  evaporation  we  would  expect  to  see  interannual  vari¬ 
ability  in  NGOC  evaporation,  which  would  have  higher 
evaporation  during  2003,  2005,  and  2007  and  lower  evapo¬ 
ration  during  2006  and  2008.  To  test  this  hypothesis  we 
analyzed  the  interannual  variability  of  the  evaporation  over 
the  NGOC  region  (Figure  6b).  Contrary  to  this  hypothesis, 
2006  exhibits  the  highest  evaporation  of  the  period  2003- 
2009.  In  addition,  there  is  a  lack  of  high-evaporation  events 
concurrent  with  the  high-salinity  events  of  2003,  2005,  and 
2007.  These  results  suggest  that  evaporation  is  not  always  the 
dominant  forcing  on  the  generation  of  the  upper  ocean  high- 
salinity  water  in  this  region  as  previously  suggested  by  Bray 
[1988],  Lavin  et  al.  [1995],  and  Lopez  [1997].  Instead, 
another  process  modifying  salinity  must  be  driving  these 


events.  In  section  3.3.2,  we  show  that  lateral  advection  of 
salinity  can  explain  these  anomalies. 

[is]  Before  moving  to  section  3.3.2,  it  is  important  to  note 
that  according  to  the  Japanese  Meteorological  Agency  El 
Nino-La  Nina  index  [Japanese  Meteorological  Agency, 
1991],  an  El  Nino  event  was  taking  place  in  the  Pacific 
Ocean  during  2006.  Concurrently,  Bray  [1988]  reported  an 
increase  in  the  NGOC  evaporation  during  El  Nino  years, 
which  agrees  with  the  high  evaporation  of 2006  of  Figure  6b. 
Although  the  evaporation  time  series  of  Figure  6b  is  in 
agreement  with  Bray’s  [1988]  results,  it  does  not  explain  the 
presence  of  the  fresh  water  during  2006. 

3.3.2.  Evaporation  and  Volume  Transport  as  the 
Generators  and  Modnlators  of  the  NGOC  Salinity 

[19]  Time  series  of  upper  ocean  transport  along  a  zonal  line 
across  the  entrance  of  the  GOC  is  shown  in  Figure  6c.  The 
2006  El  Nino  event  generated  the  second  largest  anomalous 
northward  upper  ocean  transport  over  this  period.  This 
increase  in  the  transport  is  consistent  with  the  studies  of  Strub 
and  James  [2002a,  2002b]  and  Lopez  et  al.  [2005]  who  also 
reported  an  increase  in  the  anomalous  northward  upper  ocean 
transport  into  the  GOC  during  the  1997  El  Nino  event. 
According  to  the  Japanese  Meteorological  Agency  index, 
2003  and  2005  were  classified  as  neutral  years.  However, 
those  years  include  the  strongest  upper  ocean  southward 
transport  events  of  the  period  2003-2009.  Why  does  the 
upper  ocean  northward  transport  into  the  GOC  increase 
(decrease)  during  El  Nino  (La  Nina)  years? 

[20]  During  El  Nino  (La  Nina)  years  the  equatorial  Pacific 
features  the  formation  of  downwelling  (upwelling)  Kelvin 
waves.  These  waves  propagate  eastward  until  they  reach  the 
west  coast  of  the  Americas,  where  they  excite  poleward 
propagating  coastal  trapped  waves  with  the  northward  CTWs 
reaching  the  GOC  and  beyond  [Spillane  et  al,  1987;  Enfield, 
1987;  Kessler  et  al.,  1995;  Ramp  et  ah,  1997;  Lyman  and 
Johnson,  2008].  The  along  coast  SSH  diagram  of  Figure  7 
clearly  shows  these  equatorial-originated  waves  propagat¬ 
ing  along  the  North  Americas  West  Coast  (from  the  Acapulco 
region  to  the  NGOC).  Downwelling  (upwelling)  CTWs 
induce  poleward  (equatorward)  currents  during  their  propa¬ 
gation  in  the  GOC.  Thus,  the  currents  generated  by  these 
El  Nino  (La  Nina)  induced  CTWs  strengthen  (weaken)  the 
local  PEBC  and  their  associated  upper  ocean  transport.  That 
is  displayed  by  the  along  coastal  SSS  anomaly  diagram  of 
Figure  8,  which  shows  anomalously  fresh  surface  water  being 
advected  poleward  by  the  2006  El  Nino  intensified  CTWs.  In 
concordance  with  this  explanation,  the  hydrographic  obser¬ 
vations  ofLopezetal.  [2005]  wA  Padilla  etal.  [2006]  clearly 
show  upper  ocean  low-salinity  water  in  the  NGOC  during  the 
1997  El  Nino  event  [see  Lopez  et  al,  2005,  Figure  9].  Fur¬ 
thermore,  it  is  important  to  note  that  although  the  amplitude  of 
equatorially  originated  CTWs  is  generally  increased  during 
El  Nino-La  Nina  years,  these  waves  are  also  evident  during 
neutral  years  (Figure  7).  That  is  exemplified  by  the  arrival 
of  strong  upwelling  (downwelling)  CTWs  to  the  GOC  during 
the  2003  (2008)  neutral  year  (Figure  7),  which  forced 
anomalous  upper  ocean  southward  (northward)  transport 
(Figure  6c),  and  the  consequently  maximum  (minimum)  in 
the  SSS  in  the  NGOC  (Figures  6a  and  8).  In  fact,  the  largest 
anomalous  northward  upper  ocean  transport  over  the  period 
2003-2009  occurred  during  the  neutral  year  of  2008  and  that 
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Figure  6.  (a)  Time  series  of  daily  mean  sea  surface  salinity  for  the  NGOC  region  shown  in  Figure  4. 
(b)  Time  series  of  daily  mean  evaporation  for  the  NGOC  region  shown  in  Figure  4.  (c)  Time  series  of  daily 
mean  upper  ocean  (0-100  m)  transport  (Sv)  calculated  along  the  zonal  line  across  the  entrance  of  the  GOC, 
which  is  indicated  in  Figure  Ic.  Positive  transport  indicates  northward  flow.  The  fields  were  determined  by 
an  interannually  forced  nested  GOC  HYCOM  simulation  (experiment  6  in  Table  1).  A  30  day  running  mean 
filter  has  been  applied  to  the  three  time  series.  In  accordance  with  the  Japanese  Meteorological  Agency 
El  Nino-La  Nina  index,  El  Nino,  La  Nina,  and  neutral  years  are  represented  with  red,  blue,  and  green  colors, 
respectively.  The  first  day  of  each  year  is  labeled. 
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Figure  7.  Sea  surface  height  anomaly  (color  contours  in  cm) 
time  series  determined  with  an  interannually  forced  nested 
GOC  HYCOM  simulation  (experiment  6  in  Table  1)  along 
the  west  coast  of  Mexico,  from  1 6°N  to  the  northern  Gulf  of 
California.  The  anomaly  is  relative  to  the  2003-2009  mean. 
In  accordance  with  the  Japanese  Meteorological  Agency  El 
Nino-La  Nina  index,  El  Nino,  La  Nina,  and  neutral  years  are 
represented  with  red,  blue,  and  green  colors,  respectively. 
The  first  day  of  each  year  is  labeled.  NGC  represents  the 
location  of  the  northern  Gulf  of  California. 


-0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8 


Figure  8.  Sea  surface  salinity  anomaly  (color  contours)  time 
series  computed  from  an  interannually  forced  nested  GOC 
HYCOM  simulation  (experiment  6  in  Table  1)  along  the 
west  coast  of  Mexico,  from  16°N  to  the  northern  Gulf  of 
California.  The  anomaly  is  relative  to  the  2003-2009  mean. 
In  accordance  with  the  Japanese  Meteorological  Agency  El 
Nino-La  Nina  index,  El  Nino,  La  Nina,  and  neutral  years  are 
represented  with  red,  blue,  and  green  colors,  respectively. 
The  first  day  of  each  year  is  labeled.  NGC  represents  the 
location  of  the  northern  Gulf  of  California. 
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transport  event  has  a  clear  signature  in  the  NGOC  salinity 
(Figures  6a,  6c,  and  8). 

[21]  Overall,  the  PEBC  transport  upper  ocean  low-salinity 
water  from  the  Pacific  to  the  GOC  and  that  transport  of  low- 
salinity  water  is  modulated  by  equatorially  originated  CTWs. 
Thus,  the  PEBC  are  the  provider  of  low-salinity  surface 
water.  That  water  is  used  as  the  background  conditions  by  the 
NGOC  high  evaporation  to  generate  the  NGOC  high-salinity 
surface  water. 

4.  Summary  and  Concluding  Remarks 

[22]  The  seasonal  and  interannual  variability  of  the  NGOC 
salinity  is  studied  using  an  eddy-resolving  (1/25°  horizontal 
resolution)  GOC  version  of  HYCOM  nested  in  global 
HYCOM.  The  nested  approach  (instead  of  a  pure  regional 
approach)  permits  the  free  communication  of  the  NGOC  with 
the  Pacific  Ocean.  The  NGOC  salinity  fluctuations  are 
examined  and  they  are  controlled  by  both  the  NGOC  high- 
evaporation  and  the  low-salinity  water  transported  by  the 
PEBC  from  the  Pacific  to  the  NGOC.  Model  results  indicate 
that  the  high  evaporation  generates  high-salinity  surface 
waters,  which  are  featured  by  maxima  SSS  at  the  head  of  the 
upper  NGOC  and  decreasing  SSS  toward  the  mouth  of  the 
GOC.  The  role  of  the  three  components  of  the  PEBC  (CRCC, 
MCC,  and  CTWs)  on  the  NGOC  seasonal  circulation  reversal 
is  isolated.  The  model  results  suggest  that  besides  the  wind 
and  the  annual  CTW  as  the  main  forcing  of  the  seasonal 
reversal  of  the  circulation  [Beier  and  Ripa,  1999;  Marinone, 
2003]  the  CRCC  and  the  MCC  also  contribute  to  the  gener¬ 
ation  of  the  NGOC  seasonal  circulation  reversal  (Figure  2). 
Furthermore,  the  influence  of  the  Pacific  Ocean  and  the 
Southern  GOC  (via  CRCC,  MCC,  and  annual  CTW)  over  the 
NGOC  upper  ocean  circulation  and  SSS  is  stronger  during 
spring-summer  than  during  fall-winter  (Figures  2-4). 

[23]  The  NGOC  SSS  is  characterized  by  strong  inter- 
annually  variability.  Opposite  to  our  initial  hypothesis,  the 
interannual  variability  of  the  NGOC  SSS  cannot  be  explained 
exclusively  in  terms  of  high  evaporation.  In  fact,  the  second 
lowest  SSS  of  the  period  2003-2009  occurred  during  2006, 
which  was  characterized  by  the  highest  evaporation  of  the 
period  2003-2009.  That  is  reconciled  by  the  second  largest 
upper  ocean  northward  transport  of  low-salinity  water  carried 
by  the  2006  El  Nino  intensified  PEBC  (Figures  6-8).  In  short, 
on  both  seasonal  and  interannual  time  scales,  high  evapora¬ 
tion  generates  high  salinity  and  the  PEBC  modulate  the 
salinity  by  transporting  upper  ocean  low-salinity  water  from 
the  Pacific  to  the  NGOC  (Animation  1). 
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